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Ultracold and very cold neutrons (UCN and VCN) interact strongly with nanoparticles due to the 
similarity of their wavelengths and nanoparticles sizes. We analyze the hypothesis that this inter- 
action can provide efficient cooling of neutrons by ultracold nanoparticles at certain experimental 
conditions, thus increasing the density of UCN by many orders of magnitude. The present analytical 
and numerical description of the problem is limited to the model of independent nanoparticles at 
zero temperature. Constraints of application of this model are discussed. 
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I. INTRODUCTION 

A series of experiments performed by a number of research groups have brought hght to 
the phenomenon of the quasi-elastic scattering of UCN at surfaces displaying surprisingly 
small energy changes in the order of 10"'' eV . A detailed study of this 
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process [8|, |9|, [10|, [11|, [Ij, [13| has allowed us to conclude that, for solid surfaces at least, this 
is due to the inelastic coherent scattering of UCN on nanoparticles or nanostructures weakly 
attached to the surface in a state of permanent thermal motion. This conclusion triggered 
the idea considered in this article of neutron cooling at ultracold nanoparticles. Complete 
control of the corresponding UCN losses from storage bottles and/or a significant increase in 
UCN density are of utmost importance for neutron-based research in fundamental physics. 
This research includes the measurement of the neutron lifetime Q, [3 Q| , the search 

for the non-zero neutron electric dipole moment [l^, [2^ , the study of the gravitationally 



bound quantum states of neutrons 
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23|, and the search for the non-zero neutron 
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In any experiment with trapped neutrons, nanoparticle temperature is equal to the trap 
temperature T in the typical range of 10-10^ K, while UCN energy corresponds to UCN 
temperature T mK. Ultracold neutrons therefore preferentially increase their energy in 
collisions with such "warm" nanoparticles. The probability of such inelastic UCN scattering 
on the surface is small, since the surface density of such weakly attached nanoparticles is 
typically small. 

However, the problem of neutron heating due to neutron-nanoparticle interaction can in 
principle be reversed: the interaction of warm neutrons with ultracold nanoparticles at a 
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temperature of «1 mK can cool down the neutrons [8|, |25!, |20|. If the density of weakly 
bound nanoparticles is high (these nanoparticles not only cover the surface but also fill the 
volume) and if, as the neutrons cool, the probability of their absorption and /3-decay is low, 
the neutron density will increase. This process can, for the first time, allow equilibrium 
cooling of neutrons down to the UCN temperature. One should note that in this case the 
moderator temperature should be as low as the UCN temperature in contrast to traditional 
methods to produce UCN. 

The cooling of neutrons in nuclear reactors and spallation sources by a factor of about 
10^! is achieved by just a few dozen collisions with nuclei in reactor moderators (hydrogen, 
deuterium). The energy transfer is very efficient since the mass of the moderator nuclei 
is equal to (or approximates) the neutron mass and the neutron losses during the cooling 
process are low due to small number of collisions needed to slow down neutrons. However, 
no further efficient cooling occurs: the lower the neutron energy, the larger the neutron 
wavelength. When the wavelength becomes commensurate with the distance between the 
nuclei of the moderator, the neutrons do not "see" individual nuclei any longer - they 
are just affected by the average optical potential of the medium. A further cooling of the 
neutrons due to their interaction with collective degrees of freedom (such as phonons) is 
less efficient than the moderation of the neutrons due to their collisions with nuclei. That 
is insufficient however to cool the main portion of the neutrons to the UCN energy region 
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The idea of neutron cooling on ultracold nanoparticles consists in reproducing the prin- 
ciple of neutron cooling in reactor moderators using multiple collisions. However, there is 
a difference in scale: the sizes of scattering centers are greater by a factor of w 10^; this 
increases the energy range of application of this mechanism by a few orders of magnitude. 

It should be noted that a UCN source of this type is based on the principle of UCN density 
accumulation, as in a super-thermal source [29], but not on the use of a UCN flux from 
a source in the flow-through mode. In conventional sources used to select UCN, thermal 
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equilibrium is not achieved. These sources are much hotter than UCN. Only a very small 
portion of the neutrons is used - the other neutrons are lost. Actually, these are sources of 
cold or very cold neutrons (VCN), and experimentalists have to select a narrow fraction of 
a broad energy spectrum. For instance, the most intense flux of UCN available for users is 
now produced in a liquid-deuterium source placed within the core of the high-flux reactor 
at the Institut Laue-Langevin (ILL) [28]. It increases the UCN flux by a factor of about 10^ 
in relation to that available otherwise in the reactor in the thermal equilibrium spectrum. 
Only a fraction of the neutron flux of about 10^^ is thus actually used. On the other hand, 
the cooling of neutrons on ultracold nanoparticles could provide for further neutron cooling 
in a significant energy range, thereby increasing the neutron density. 

The new method for producing UCN consists in the equilibrium cooling of VCN - through 
their many collisions with ultracold nanoparticles made from low-absorption materials (D2 , 
D2O, O2 etc.) - down to the temperature of these nanoparticles of « 1 mK, during the 
diffusion motion of these neutrons in a macroscopically large body of nanoparticles. The 
principle of equilibrium cooling allows an increase in the neutron phase-space density, in 
contrast to the method of selecting a narrow energy range out of a warmer neutron spectrum. 
The use of nanoparticles provides a sufflciently large cross section for coherent interaction 
and an inhomogeneity of the moderator density, on a spatial scale of about the neutron 
wavelength; it also shifts the energy transfer range far below a value of about 10^"^ eV, 
the characteristic limit for liquid and solid moderators. Many collisions are needed since 
the mass of the nanoparticles is much larger than the neutron mass; the energy transfer to 
nanoparticles and nanostructures is only moderately efflcient. The need for a large number 
of collisions limits the choice of materials: only low absorption materials are appropriate. 
The temperature of the nanoparticles must correspond to the minimal energy to which 
neutrons can still be cooled using this method. The diffusion motion of neutrons in the 
body of nanoparticles allows us to minimize the thermalization length and, accordingly, to 
increase the achievable UCN density. 

The cooling itself is provided by the interaction of neutrons with individual degrees of 
freedom of weakly bound or free nanoparticles, as well as by the excitation of collective 
degrees of freedom in the body of nanoparticles (e.g. vibrations and rotations), and also 
by the breakin g of inter-particle bonds. Details about gels of nanoparticles can be found 
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36|. Even free nanoparticles in the gel have several degrees of freedom: 



rotation and translation. In this paper, we provide detailed calculations of the cooling 
of neutrons in a gel of nanoparticles, considering only the collisions on free nanoparticles, 
and neglecting the rotation. The interaction between the nanoparticles, including the long- 
range interaction induced by helium, is neglected here. So we deal with an idealized gas of 
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free nanoparticles at K in superfluid helium. Our goal is to determine the behavior of a 
neutron in such a moderator. 

In section|lIl we give a model for the interaction of a neutron with a single nanoparticlc in 
suspension in liquid helium. This model leads to a complete quantum solution for three main 
quantities, which are, the absorption cross section (Ja, the total scattering cross section as, 
and at last, the mean relative energy loss per collision ^ = This allows us to describe 

the slowing down of neutrons in a gas of free nanoparticles. In section IIIH we deal with 
the ideal situation of the infinite moderator, an infinite medium made with nanoparticles 
surrounded by helium at K, and where the only loss of neutrons are due to the absorption 
by a nanoparticlc (/3-decay is neglected). This naive model provides necessary conditions 
for an efficient moderation. In section HVl we will present estimations for finite moderators, 
the characteristic size and the characteristic time of thermalization. We compare efficiency 
of the moderator with nanoparticles made of different materials, and also for different 
nanoparticle sizes. We first chose deuterium as the material to illustrate the calculations. 



II. FREE NANOPARTICLES MODEL 



Let us consider a nanoparticle with radius R, made of hundreds of nuclei, immersed into 
superfluid helium. A low energy incident neutron only sees the average potential of each 
nanoparticle. So we will assume the following phenomenological potential for the interaction 
of a neutron with the considered nanoparticle: 

( V = Vo-iVi ifr<R, 
Vir)=\ (1) 
[ if r > i?. 

We will provide detailed calculations of the collision parameters in the Born approximation 
using this model. But let us first estimate all the parameters of the potential ([T]) describing 
neutron-nanoparticle interaction. 



Parameters of neutron-nanoparticle interaction potential 

The potential Vq — iVi is taken to be the averaged of each nucleus Fermi potential in the 
nanoparticle (neutron- helium interaction potential has to be subtracted): 

Vo -iVi = Yl P^'^ J (^)^^ - ■^H- (2) 
j 

where the sum is done over the different kinds of nuclei in the nanoparticlc - for instance, for 
a heavy water nanoparticle, we must take into account the contributions of both deuterium 
and oxygen nuclei. In each term of the sum, p^^^ is the number of nuclei of type j per 



5 



TABLE I: Nuclear data used in these calculations. The experimental data are taken from 



Nucleus b (fm) aa (mbarn) Mass (m) 

-3.74 333 0.999 

6.67 0.52 1.997 

^He 5.33-10^ 3.968 

■*He 3.26 3.968 

Be 7.79 7.6 8.935 

C 6.65 3.5 11.91 

O 5.80 0.19 15.86 



volume unit, and V^^^ = Vq — iV} is the interaction potential between a neutron and a 
j-type nucleus. The real part of this potential can be derived from the coherent scattering 
length fo^J) 

where m is the neutron mass. The imaginary part of the potential, which describes the 
possibility of the neutron capture, can be calculated from the absorption cross section via 
the optical theorem, taking into account the fact that this interaction is not strong and can 
be treated within the first Born approximation: 

|J Jvl^Ur^a^^Hk)k. (4) 

The experimental data are given for thermal neutrons, i.e. for a neutron velocity of 2200 
m/s. So we express the result in term of (Ta^ (A;o)fcoj where fcg is the wave vector of a thermal 
neutron. Eventually: 

Vo = 27r-Vp«6«, (5) 
m ^-^ 

3 

Vi - |^^p«aW(fco)fco. (6) 
Zm ^ — ' 

i 

From the nuclear data shown in table |T1 we can compute all properties of the potential of 
interaction between neutrons and nanoparticles in superfluid helium. The results are shown 
in table ini for various materials, taking into account the following points: 

• We note A the nanoparticle mass in units of neutron mass. Its variation with the 
radius is determined by the parameter Aq: A = Ai^R^ . 

• The density of liquid helium is 124.9 kg m~'^ at boiling point 4.25 K. We calculate 
the effective potential of the nanoparticles in liquid helium which has a potential 
Fhc = 15.9 neV. 
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TABLE II: Nanoparticle characteristics. 



Nanoparticle 


Density 


M 


No 


Vo 


Vi 




(kg/m^) 


(nm-^) 




(neV) 


(neV) 


D2 


195 


488 


0.001510 


85 


2.2 -10"^ 


D2O 


1020 


2551 


0.000292 


137 


2.7 -10"^ 


O2 


1230 


3076 


0.000242 


54 


0.64 •10"'^ 


CO2 


1560 


3901 


0.000191 


85.5 


6.0 -lO-*^ 


C (Diamond) 


3520 


8803 


0.000085 


290 


45 •10"'^ 


Be 


1850 


4627 


0.000161 


235 


68 •10"'^ 



• We calculate the density of nanoparticles (number of nanoparticles in a unit volume) 
N — NqR^^, assuming that the total mass of the nanoparticles is 1% of the total mass 
of the helium. This parameter is only useful for the description of a finite moderator 
in section ITVl 

Validity of the model at high energy 

At high velocities neutrons can see individual nuclei, and the form of the potential we 
assumed is not valid. This happens when the wavelength of the neutron is smaller than the 
inter-atomic distances. So we use the following limit for our model: 

27r , , 

Amin = T = « (7) 

"'max 

where d is the mean inter-atomic distances, and fcmax is the maximum neutron wave vector 
allowed by our phenomenological model. The table shows this limit for various materials. 



Nanoparticle 


D2 


D2O 


O2 


CO2 


C Be 


d (nm) 


0.32 


0.32 


0.35 


0.36 


0.18 0.20 


fcmax (nm^i) 


19 


19 


16 


17 


35 31 



As we will see, this constraint is not of actual importance in the range of parameters 
concerned, corresponding to the optimal conditions of neutron thermalization. 

There is another limit at high velocities. The neutron can excite the internal degrees 
of freedom of nanoparticles, so that the collision can be inelastic. We expect that the 
probability of phonon excitation is low because the phonon wavelenght should be in this 
case shorter than the considered size of nanoparticles. Also this process can only increase 
the efficiency of the neutron cooling. 
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Scattering amplitude within Born approximation 

Our goal is to calculate the three relevant quantities describing the collision of a neutron 
with a nanoparticle: the total scattering cross section, the absorption cross section, and the 
mean energy loss. Let us note that in this problem the absorption probability is extremely 
small with respect to the elastic scattering one. Therefore the total interaction cross section 
is approximately equal to the elastic one. 

In the first Born approximation we can easily compute this scattering amplitude in the 
center-of-mass system (c.m.s.), that is, the amplitude for a neutron with incident wave 
vector k to be scattered at final wave vector k': 

fie) = -^^-^ [e^^^^^'y^V{r)dr (8) 



An h? _ 

being the angle between k and k'. Let q = k — k' be the momentum transfer. The 
collision is elastic in the c.m.s. so that k = k' and 

g = 2fcsin(6i/2). (9) 

Finally, we find 

9m 1 

where Ji(^) is the first spherical Bessel function: 

sin(X) cos{X) 



= (11) 



^From the scattering amplitude we can calculate the elastic cross section: 
as^ j |/|'dr! = 2^ 



2rn 



2kR ^ 



I{kR) = / -ji{xYdx 







X 



\ ( ^ 1 sin(4A:i?) sin^(2fci?)\ 



where 



4 y (2kRY ^ [IkRf (2kRY ) ■ ^^^^ 
By using the optical theorem, we can calculate the absorption cross section aa at the first 
order Born approximation: 

lm{f{e = 0)) = (14) 
As /(O) = ~l^VR^, one obtains: 

47r 2m „4 1 , ^ 

- T 7^^^ ^ kR- (1^) 

Calculated elastic and absorption cross sections are presented in Fig. [T] as a function of 
neutron velocity. The calculations are performed for different nanoparticle's radii. For low 
energies where kR ^ 1, one can easily see from this figure as well as from (|T2|) and ^T5\i 
that as R^ and (Ja ^ R^- 
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FIG. 1: Elastic and absorption cross sections as a function of neutron velocity, for three values of 
the deuterium nanoparticles radii; 1, 2, and 5 nm. 



Mean energy loss per collision 



We consider a collision between the neutron and a nanoparticle which has a mass A, in 
units of the neutron mass. Let 9 be the scattering angle. As the mass of the nanoparticle is 
much greater than the mass of the neutron, we know that the energy transfer, in the first 
order in l/A is given by: 



AE 



~4lsin(0/2)2 



(16) 



Since we know the scattering amplitude at the Born approximation, we can calculate the 
mean relative energy loss: 



{AE) _ J^§^dn _2n r AE 
E I^dn aJo E 

In terms of the variable x = qR = 2kRsm{0 /2): 

AE 



sin 



E 



1 2m ,1 



Now we can express the mean relative energy loss: 

(AE) _ 1 1 J{kR) 



E 



A [kRY I{kR) 



with 



J{kR) _ Q^"' xji{xYdx 



(17) 



(18) 



(19) 



(20) 
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FIG. 2: Relative mean energy loss per collision ^, as a function of neutron velocity, for three values 
of the deuterium nanoparticles radii: 1, 2, and 5 nm. 



The first integral in this expression can be expressed in terms of special functions but has 
to be calculated numerically where the second one was calculated previously in (|13p . 

The relative mean energy loss as a function of neutron velocity is presented in Fig. [2] for 
different values of nanoparticles radii. Let us note that, for small energies where kR ^ 1, 
this loss behaves as ^ ^ R^^ . 



Interpretation of the results 

As expected, the absorption cross section follow a l/w law, and is proportional to the 
nanoparticle's mass vloi?"^. At low velocities, the elastic cross section and the energy loss 
are constant up to some critical momentum kc ~ beyond which the collision becomes 
anisotropic. The total elastic cross section and the energy loss therefore decrease rapidly. 

As we have mentioned previously, the probability of collision is strongly dependent on 
the size of the nanoparticle, as R^: due to the coherent scattering of a neutron at nuclei in 
the nanoparticle, the scattering cross section is proportional of the square of the number of 
nuclei in the nanoparticle. From the condition kc ~ the critical momentum is higher 
for smaller nanoparticles. 

The efficiency of the moderation results from competition between the absorption process 
and the elastic scattering. The probability of the first process is proportional to the absorp- 
tion cross section a a multiplied by the number of collisions needed to slow down neutrons 
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(this number is inversely proportional to the averaged energy losses ^). The probability of 
elastic scattering is proportional to the elastic cross section a^. As we have emphasized, all 
these three values da, S,, <^s are strongly i? dependent but their combination Gaji^Os which 
defines the moderator efficiency is not. This combination depends only on the product fci? 
and has its minimum for fci? ~ 1 (as it can be seen from Fig. [T]and Fig. [2]). Therefore it is 
natural to expect the best moderator properties for neutron velocities close to the critical 
one 1/R- This ratio decreases for both higher and lower velocities. We will see that a 
more detailed analysis in section Hill will confirm this general statement. 



Validity of the first Born approximation 

Sufficient condition for the validity of the first Born approximation is: 

'^-^V.R' « 1 (21) 

The lower the radius, the better the approximation. To see where this condition is ver- 
ified, let us calculate the maximum radius i?Born to have an accuracy better than 10%: 



Nanoparticle D2 D2O O2 CO2 C Be 



i?Bom (nm) 4.5 3.7 5.4 4.5 2.6 2.9 

A 



Actually, the approximation works much better, because our potential is repulsive 
more accurate calculation, with partial waves given in Appendix [A] shows that the accuracy 
of the Born approximation is sufficient enough for our needs. 



III. MODEL OF THE INFINITE MODERATOR 



The role of any moderator is to increase the density of neutrons in the phase space. In this 
section, we will determine the influence of the properties of nanoparticles - their material 
size - on the evolution of neutron density. In order to achieve that, we have to identify the 
criteria for cooling efficiency. 

The most general problem that we perceive is the following. We have a moderator medium 
made of free nanoparticles in suspension in supcrfluid helium. There is a neutron source, at 
a certain location, with a certain velocity distribution. What is the neutron density at every 
point in the moderator medium and what is the distribution of velocities? The quantity 
describing the state of the system in this problem is a seven-variable function, the density 
of the neutrons in the phase space n(r, k, t). In general case, this density evolves according 
to the transport equation, or Boltzmann equation, and all the variables can be coupled. 
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Solution of this equation in general case can be quite complicated task so we will propose 
a first necessary criterium of cooling efficiency within the model of infinite moderator. In 
this model, we eliminate the space variable in the density; we assume that the sources are 
uniformly distributed, and that the medium is homogeneous and isotropic. The quantity 
which describes the system is only a function of the energy (and also of time, but it is only 
the stationary regime that interests us) . We can analytically compute the energy spectrum 
of the neutrons in this moderator, using the three quantities Cs, a a, C = "^^r^ describing 
the elementary processes. 



The energy spectrum in the infinite moderator 



To compute the neutron energy spectrum in the moderator n[E,t), we have to use the 
equation for the conservation of the number of neutrons: 

— - ^ + Absorption - Source = (22) 

where q{E, t) is the cooling current, i.e. the number of neutrons scattered from an energy 
greater than E to an energy lower than _B, in the unit of time. The absorption term 
can be related to the macroscopic absorption cross section = Nua via the flux variable 
4>{E, t) — vn{E, t) = ^J^n{E, t), so that the conservation equation, using the flux variable, 
can be expressed as: 

m dcf) dq 



^ . Eacj) = Source (23) 
2E dt dE ^ ^ ' 

Now we have to specify the cooling current under the assumption that ^ 1 which is clear 
to be a very good approximation for neutron scattering on nanoparticl es ( see Fig. [5]) . Under 
this assumption the cooling current has the simple and natural form 39|: 



q{E) = EaE)J:s{E)^{E). (24) 

Let us consider the stationary regime. The conservation equation is completely expressed 
in terms of our three microscopic quantities, and so the time independent equation is given 

by 

- ^ [^^sE(P] + Sa<A = (25) 

We have not expressed the source term, because if the source is punctual (a given flux <I>o 
at a given energy Eq) then we can solve the equation without any source in the domain 
< E < Eq and the source is a boundary condition: 4i{Eq) — ^q. So now, as we have a 
usual first order linear differential equation, we can solve it analytically with the previous 
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FIG. 3: Phase space density in stationary regime as a function of velocity, for three radii of deu- 
terium nanoparticles. The source is assumed to be monochromatic, with initial energy equal/greater 
than the maximum energy presented in the graph. 



boundary condition: 



We can notice that, if there is no absorption, and if both ^ and Eg are energy independent, 
we find the usual behavior, the so called the 1/E flux law. 

From the flux, it is useful to derive the density n(k) of neutrons in the phase space; 

n{k) = n{ko) — exp -2 / — (27) 



kj (es,)(fc) " \ A- « 



A necessary condition for efficient cooling 



The density of neutrons in the velocity space, for the stationary regime in the infinite 
moderator, is presented in Fig. [3l 

To interpret this figure, we have to assume that there is a monochromatic neutron source 
at velocity fco, and we extract the neutrons at the velocity fci < ko. Then the relative 
increase of density in the phase space is given by n{ki)/n{ko). As the goal is to increase the 
density of neutrons in the phase space, we see that the curve gives us a necessary condition 
for the efficiency of the cooling. Indeed, the cooling is efficient only if n(fci) > n{ko), that 
is in the decreasing part of the curve. This differential necessary condition is then given by 
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FIG. 4: Efficiency domain of the moderation as a function of nanoparticle's radii, for different 
materials. Interference effects are not taken into account. 

^^p(fco) < 0, that we can put into the form: 

E„ ldln(^ 

The coohng domain is the velocity domain in which this condition is satisfied. Figure 3] 
shows that the cooling is not efficient both at low velocities and at high velocities. 

At low velocities, the cooling is not efficient because the absorption is important com- 
pared to the diffusion cross section and the energy loss, and at high velocities (basically at 
velocities higher than the critical velocity), the cooling is not efficient because both the scat- 
tering cross section and the energy loss fall down. Figure 2] shows the cooling domain as a 
function of the radius of the nanoparticles, for several materials. These curves leads to three 
remarks. For certain materials, such as diamond, and for certain radius of nanoparticles, 
there is no cooling domain at all. The cooling domain is bigger for smaller nanoparticles, 
because for smaller nanoparticles, the critical velocity is higher. The domain for deuterium 
nanoparticles include the domain for all the other materials. 

We can now put a limit for the total increase of density in the phase space achievable for 
the infinite moderator and for a monochromatic source. Indeed, in Fig. [3] we see that the 
best we can do is to have a source with fco the minimum of the curve, and to extract 
neutrons with fci at the maximum of the curve. The maximum gain is given by n{ki)/n{ko). 
Figure [5] shows this limit, as a function of the radius of the nanoparticles, and for different 
materials. We see that this maximum is bigger for smaller nanoparticles, because the cooling 
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FIG. 5: Maximum compression of phase space density possible, as a function of the radius of the 
nanoparticles, and for different materials. Interference effects are not taken into account. 

domain is higher. The best material for this criterium is deuterium. 

IV. MORE REALISTIC MODERATORS 

The consideration of the infinite moderator gives us a necessary condition for the efficiency 
of the cooling. The results presented do not depend on the density of nanoparticles in the 
moderator. When it comes to realistic moderators, we have to take two characteristics of 
the moderation into account: the size of the moderator, and the thermalization time. These 
two quantities depend on the density of nanoparticles, and for the calculations we assumed 
that in a given volume, the total mass of nanoparticles represents 1 % of the total mass 
of helium. That means that for a gel of smaller nanoparticles, the distance between the 
nanoparticles is assumed to be smaller. 

Let us pursue the idea to estimating the size of the moderator and the thermalization 
time. Suppose a monochromatic neutron source at velocity fco inside the efficient domain, 
and that we want to increase the density in phase space by the factor e. According to 
Fig. [3l we have to extract neutrons in the moderator with an energy fco — Afc, where Afc 
corresponds to an increase of factor e for the function n(fc): 

(29) 



Afc = 



dlnr7,(fc) 
dfc 



We know that this decrease of velocity by Afc corresponds to a certain number of collisions 
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AN 

AN^l^ (30) 

During AN collisions, neutron assume a Brownian-like trajectory, and we define L^k^) as 
the square-mean-root of the distance travelled. We also define the thermalization time T(fco) 
as the time it takes for the neutron to collide AN times. 

Estimation of the moderator's size 

The square mean distance travelled by a neutron Ar^ after A collisions is given by [4^ ■ 

Ar^ - J^AN (31) 

This expression was obtained under the assumption that the diffusion cross section and 
the mean relative energy loss £ are energy-independent, and that the scattering is nearly 
isotropic. 

From equations ([29]) and ((30|) . we can estimate the characteristic length L{k) of the 
moderator, corresponding to an increase of factor e of the phase space density: 

1 dln(n) A^^^ 2dln(n) 



L{kY dr2 4''- " dln(fc) 

ldln(eS,) 



1 



(32) 



2 dln(fc) iT.s 

Note that we find our necessary condition ([28]) . in the sense that the condition is satisfied 
if and only if the size L is real. 

Estimation of the thermalization time 

Now we must estimate of the thermalization time. We know that the mean time between 
two collisions is given by: 

h kljs 

The time needed for a neutron to collide AA^ times is therefore simply TcouAiV. From 
equations (j29p and (j30p . we can estimate the thermalization time r(fc), corresponding to an 
increase by a factor e of the phase space density: 

1 ^ j2t^ dln(n) 



T[k) 2m dfc 



hk 

2m 



^^ldln(eE,) Ea 



(34) 



2 dln(fc) 

Again, we find our necessary condition (l28l) . in the sense that the condition is satisfied if 
and only if the characteristic time r is positive. 



16 




0.01 ' ' ' ' ' ' ' 

1 2 3 4 5 6 

R (nm) 

FIG. 6; Moderator size needed to increase the density in phase space by a factor e. Interference 
effects are not taken into account, (a) Moderator size as a function of velocity, for three radii 
of deuterium nanoparticles. (b) The minimum moderator size as a function of the nanoparticles 
radius, for different materials. 



Comments on the results 



The results are shown in Fig. [S] and Fig. [71 Let us first comment on the figures for 
the deuterium example. We see first that outside of the efficiency domain of the cooling, 
both the thermalization time and the moderator size diverge. We also remark that the 
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FIG. 7; Thermalization time needed to increase the density in phase space by a factor e. Interference 
effects are not taken into account, (a) Thermalization time as a function of velocity, for three radii of 
deuterium nanoparticles. (b) The minimum thermalization time as a function of the nanoparticles 
radius, for different materials. 



moderator size is nearly constant over the efficiency domain, although the thermalization 
time presents a clear minimum near the critical velocity. But the main difference between the 
two conditions is the dependence on the radius of the nanoparticles. For small nanoparticles, 
we need fewer collisions to cool down the neutrons, and wc can see in Fig. [7] that the 
smaller the nanoparticles, the smaller the thermalization time. On the other hand, for 
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small nanoparticles, the mean free path is high, and we see in Fig. [S] that the smaller the 
nanoparticles, the bigger the moderator size. We cannot therefore optimize the radius of 
the nanoparticlc in a general way; we must look at more practical considerations. 

We can use both the size of the moderator and the thermalization time as criteria to 
compare materials. We can see that using both criteria, deuterium is better than any other 
material by almost one order of magnitude. We also see that the absolute values are very 
competitive. Indeed, to increase the density in phase space by factor e, the characteristic 
size is a few tens of centimeters and the thermalization time is less than one second, for a 
typical radius of 2 nm. 

Practical proposal 

To achieve a compression in the phase space by many orders of magnitude, we can as- 
sociate the nanoparticle's moderator to a compression in real space. To explain this idea 
in more details. Fig. [8] may be of help. We start with a gas of neutrons uniformly located 
in the initial volume Vq with a velocity distribution with a mean value vq and a width Av. 
We decrease the volume by a factor V/Vq. According to the Liouvillc theorem, the width 
of the velocity distribution will increase as Av — > (^)^/^Aw, so that the density in the 
velocity space decreases by the same factor (in absence of absorption this relation would 
be precise, it is approximately valid in the domain of efficient cooling of neutrons defined 
above). This is the case when there is no moderator, but if we put a nanoparticle's modera- 
tor in the volume (not in all the volume, but in a little box in a corner), the moderation will 
compensate the increase of density in the velocity space. So finally the maximum density 
in the phase space achievable is the product of the initial density, the volume compression 
Vo/V, and the maximal compression due to the moderator shown in Fig. [5l As we can in 
principle reach as much volume compression as we wish, there is no theoretical limitation 
of the maximum compression in phase space. With this method, it is not necessary to use 
all the cooling domain shown in Fig. [5] (two orders of magnitude of velocity in a favorable 
case); we can concentrate the process around the optimal velocity. This optimal velocity 
should be chosen according to two practical considerations: 

1. We are limited by the size of the moderator, which should not exceed a few tens of 
centimeters for practical purposes. 

2. The optimal velocity should correspond to a temperature greater than that achievable 
with a dilution cryostat, which is about 10 mK. The corresponding wave vector is 
about fco = 0.2 nm~'^. 
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V 

10 mK 

t = Nanoparticle's moderator t = T 

FIG. 8; Practical proposal to increase the neutron density in phase space, combining a compression 
in the real space with a piston as well as a compression in the velocity space with the nanoparticle's 
moderator. 

This two requirements are in competition, so we propose a practical criterium defined as : 

A )l3 1,3 _ l3 

-JT- - Jl • 

It is the ratio between the phase space volume for which we can have an efficient cooling in 
an actual device and the characteristic volume of the moderator. This criterium is plotted 
in Fig. [9] for deuterium nanoparticles as a function of the size of the nanoparticle. This 
figure shows that the interresting range for nanoparticles radius is between two and five 
nanometers. Althought the thermalization time is of principle importance, we showed that 
for deuterium nanoparticle this time is in any relevant cases small compared to the neutron's 
lifetime. 




V. EXTENSIONS OF THE MODEL 



For this special case of moderator nanoparticle gels there are several limitations to the 
validity of our model. One obvious limitation is that chains of nanoparticles are not taken 
into account; in this paper we neglect such chains and consider only the cooling of neutrons 
caused by the collisions with free nanoparticles. We have already mentioned the limitation 
at high energies in section this limitation is inherent to our assumptions about the 
interaction of a neutron and a single nanoparticle. But there is a limitation also at low 
velocities, which is specific to the gel. 
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Limitation at low energy 



At low velocities, we cannot neglect the interferences between the waves diffracted by 
several nanoparticles. This happens when the wavelength of the neutron is bigger than the 
distance between the nanoparticles. We see that this limitation does not come from our 
assumptions about the interaction of a neutron with a single nanoparticle. This limitation 
is specific to the practical medium that we plan to use, that is, the gel of nanoparticles. 
Let D be the distance between the nanoparticles in the gel; we can compute D, assuming 

1 /3 

that the total mass of nanoparticles is 1% of the total mass of helium: D = R/N^ . In 
the following table we give D for nanoparticles with radius 1 nm. We also give the minimal 
wave vector for the validity of the model k^in, defined as follow: Amax — -jr^ — D. 



Nanoparticle D2 


D2O 


O2 


CO2 


C Be 


Di am (nm) 8.7 


15.1 


16.0 


17.4 


22.7 18.4 


fcmin (nm^i) 0.7 


0.4 


0.4 


0.4 


0.3 0.4 



We can see that this is a serious limitation, but not that dramatic. Actually, for prac- 
tical purposes, the lowest energy we are interested in is that corresponding to the lowest 
temperature achievable with a dilution cryostat, which is about 10 mK. The corresponding 
wave vector is about 0.2 nm^^. So we only need to know , for practical purposes, the 
first correction of , <^a and ^ due to the interference of neutron waves on the neighboring 
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nanoparticles. And this can be done, because the first correction can be estimated consid- 
ering only two nanoparticles. This approach is valid until the neutron wavelength covers 
three or more nanoparticles, that means that it is valid until fcniin/2. 

Let D be the vector distance between two nanoparticles. At the first order of Born 
approximation - neglecting multiple diffusions - the scattering amplitude for a neutron 
colliding on this system of two nanoparticles is given by: 

/(D,q)^(l + e*'i°)/(q) (36) 

where q is the momentum transfer and /(q) is the scattering amplitude for the collision on 
a single nanoparticle. From that result, we can conclude that the amplitude for the forward 
scattering - q = - is simply twice the amplitude calculated for a single nanoparticle. We 
can conclude that the absorption cross section is not affected by the interferences. Now, to 
calculate the effects of interferences on the scattering cross section and the energy loss, the 
physical relevant quantity is the average differential scattering cross section, averaging on 
all possible directions for D: 



= /l/(D,q)| 



|1 + „*0|^^D/D , 

47r 



s'm{qD) 



2(^l + ^^^^jl,/(q)l' (37) 

Is is now possible to calculate as and using this new effective differential cross section. 
Figure [To] shows our plot of the relevant quantity, ^CTs, we can call the energy loss cross 
section, and the one calculated with the estimated correction. One can see that the first 
order correction is not important. 



Purity of Deuterium and Helium 

In all our calculations, we assumed that the helium was purely ^He, with no absorption, 
therefore, by the Helium medium. We also assumed that the deuterium in the nanoparticles 
was pure. Let us estimate the purity we actually need. 

Firstly, consider the presence of hydrogen ~ which is an efficient neutron absorber - inside 
the deuterium. We can compute the purity needed so that the absorption differs from the 
ideal case by less than 10 %. If x is the proportion of hydrogen relatively to deuterium, 
then the absorption cross section is: 

CTaix) = (1 - x)aa{x = 0) + Xffaix = 1). (38) 

Using the table U we can find that the requirement of less than 10 % increase is satisfied if 
a; = 1.5 • 10~^. We can easily achieve this purity. 
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FIG. 10: is plotted as a lunction of velocity, for a deuterium nanoparticle of radius 1 nm. The 
dot line is calculated without interference effects, and the continuous line take into account the 
first correction due to interference effects. 

Let us now consider the presence of "^He. We can estimate the hfetime Ta of a neutron due 
to the absorption by '^He. If x is the proportion of ^He, (T^ is the absorption cross section for 
thermal neutrons (with velocity vq = 2200 m/s), and A^hc is the density of hchum nuclei, 
then: 



If we have a purity x = 0.5 • 10^^ , the absorption hfetime Tq is ten times the intrinsic 
neutron hfetime. The natural abundance of '^He is 1.4 ■ 10^^, and the purity we need is 
achievable. 



A new concept for producing high UCN density is analyzed within the framework of 
the free nanoparticles model. This concept is based on neutron cooling using ultracold 
nanoparticles of deuterium, heavy water, etc. We have shown that increase in the phase 
space density of neutrons, within the model of free nanoparticles, is possible, given certain 
parameters of nanoparticles and neutron velocity. Thus, solid deuterium, which is shown 
to be the best material, provides efficient cooling of neutrons in the range 1 - 10^ nr/s, in 
an infinite medium of free nanoparticles of radius 1-2 nm, sufficiently spatially separated. 




X NucCTaVO 



1 



(39) 



VI. CONCLUSIONS 
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The characteristic cooling time is much shorter than the corresponding absorption time, or 
the neutron /3-decay Hfetime for optimum parameters of neutrons and nanoparticles. The 
moderator size of, at most, a few times 10 cm ahows in principle the realization of the 
cooling mechanism presented. 

We examine the different constraints on the model, such as scattering at individual nuclei 
(for too short wavelengths) , excitation of internal degrees of freedom in a nanoparticle (such 
as phonons), extensions of the Born approximation description used (the partial waves ex- 
pansion), neutron optical effects due to diffraction at several nanoparticles simultaneously, 
the purity of deuterium or helium. We show that these constraints do not change our main 
conclusions. We do not consider in this present article such effects as rotation of the nanopar- 
ticles, interaction between nanoparticles, in particular excitation of the collective degrees of 
freedom for nanoparticles in gels or any breaking of the intcr-nanoparticlcs bounds. We do 
not consider neither the influence of the non-zero temperature of the nanoparticles on the 
cooling process. These phenomena are expected to be analyzed in further publications. 



APPENDIX A: PARTIAL WAVE EXPANSION 



V (m/s) 




0.01 0.1 1 

k(nm"'') 

FIG. 11: Elastic and absorption cross-sections as a function of neutron velocity, for a deuterium 
nanoparticle of radius 5 nm. The two bold lines are calculated using the first Born approximation, 
and the thin lines are calculated using the partial wave method, till the order 0, 1, 5 and 9. 

Since the potential interaction between a neutron and a nanoparticle is approximated by 
spherical rectangular barrier, it is possible to give exact solution for the scattering amplitude. 
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using the partial wave expansion. 

The analytical solution for the scattering amplitude rji for /th partial wave is known 
to be equal to 



7; h-{kR) - k h- (kR)) 
k hf'{kR)--ii h+{kR) 



where hf are spherical Hankel functions and 



ji{KR) 



where ji is spherical Bessel function and — k"^ — 2to(Vo + iVi)/fi^ . 
The elastic and absorption cross sections are then given by 



^ ^E(2' + l)|l-^'l'' (A3) 

oo 
1=0 

The result is shown in Fig. 111! for a deuterium nanoparticle of radius 5 nm. We first 
estimated this size to be the limit of application of the Born approximation. The figure 
shows that the Born approximation is actually precise enough. The width of the bold line 
in the figure corresponds to a relative error of 5%, and there is no deviation between the 
born approximation result and the partial wave result at this precision level. 
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